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ABSTRACT 

We present detections of numerous 10-20/im H2O emission lines from two 
protoplanetary disks around the T Tauri stars AS 205A and DR Tau, obtained 
using the InfraRed Spectrograph on the Spitzer Space Telescope. Follow-up 3- 
5 fim Keck-NIRSPEC data confirm the presence of abundant water and spectrally 
resolve the lines. We also detect the P4.5 (2.934 /xm) and P9.5 (3.179 /im) dou- 
blets of OH and ^^CO/^^CO f=1^0 emission in both sources. Line shapes and 
LTE models suggest that the emission from all three molecules originates between 
~0.5 and 5 AU, and so will provide a new window for understanding the chemical 
environment during terrestrial planet formation. LTE models also imply signifi- 
cant columns of H2O and OH in the inner disk atmospheres, suggesting physical 
transport of volatile ices either vertically or radially; while the significant radial 
extent of the emission stresses the importance of a more complete understanding 
of non-thermal excitation processes. 
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Introduction 



One of the most intriguing questions in the study of the formation of planets, and of 
terrestrial planets in particular, is how water is transported to their surfaces, and whether 
or not water is a common ingredient during their formation and early evolution. Spatially 
and spectrally resolved observations of water in extrasolar planetary systems and protoplan- 
etary disks would be instrumental in resolving these q uestions. Unfortunately , detecting 



water in extrasolar planetary systems remains difficult (lEhrenreich et aL 



in protoplanetary disks has be en elusive, with o nly a few definitive ice (IMalfait et al.lll998 



200 



3) and 



water 



Terada et al.l 120071 ) and vapor (ICarr et al.ll2004j ) detections. Recently, however, the Spitzer 
InfraRed Spectrograph (IRS) has begun to reveal warm water vapor emission lines from 
disks around classical T Tauri stars (cTTs, Carr & Najita 2008). 

In this Letter, we present high-resolution spectroscopy of numerous water vapor emission 
lines from protoplanetary disks around two cTTs: P R Tau and AS 205 A. Spectra were 
obtained using both the Spitzer-IRS and NIRSPEC jMcLean et all llQQsj ) at the Keck II 
telescope. Accretion rates onto these stars are high, which causes the inner disks to be 
significantly hotter than in more quiescent systems. DR Tau is known to have variable mass 
accretion, with measured ra tes between 0.3 and 79 x 10^"^ MqJt^^ (iGullbring et al.l l2000l : 
Johns-KruU fc Gaffordll2002h . AS 205A has an accretion rate of 7.2 x lO^^MQyr"^ and is 
the prir nary component o f a 1.3" triple system, with the secondary being a spectroscopic 
binary (lEisner et al.l l2005l ) . However, it is likely that the physical separation of AS 205 A 
and B is significantly larger than the pr ojected separation and that , therefore, AS 205B does 
not affect the (inner) disk of AS 205A (lAndrews fc WiUiamsl 120071 ). 



Observations 



The Spitzer IRS Short-High (SH) r nodule spectrum of AS 205A was taken as par t of the 



"cores to disks" (c2d) Legacy program (lEvans et al.ll2003l : iKessler-Silacci et al. 



2006), while 



the sp ectrum of DR Tau, observed as part of the guaranteed time observations (IHouck et al. 
2004 ). was obtained from the Spitzer archive. The 2D ba s ic cal ibrated data images were 
reduced using the c2d pipehne as described in iLahuis et al.l (120061 ). 



To confirm the detections as well as measure the gas emission line profiles, we obtained 
high resolution (A/AA=25000) L-band spectra using NIRSPEC in a single grating setting, 
of which two orders, centered at 2.92 /im and 3.16/im, respectively, have identifiable H2O 
and OH emission features. We also collected CO t'=l — »0 M-band spectra (~5/im) as part 



of an on-going survey (c.f. iBlake &: Boogertll2004l . ISalyk et al.ll2007l ). The data were reduced 
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in a standard way. Division by the spectra of standard stars (HR 1620 (A7V) for DR Tau; 
HR 5993 (BIV) for AS 205 A) was used to correct for telluric absorption, which can be 
poor when the absorption is high. Thus, regions of the spectrum with <70% atmospheric 
transmission were removed. Fluxes were calibrated using the photometry and spectral types 
of the standards, and are generally accurate to 



0. 



The NIRSPEC line s are narrow (FWHM5, 35 km s~^) , spatially unresolved, and cen- 
tered at the stellar ulsr ( iHerbig fc Belllll988l : lThi et al.ll200ll ) to within our uncertainty (a few 
km s~^) — all consistent with disk emission. Because the disks are optically thic k and heated 
from above, the line emission likely arises from the upper disk atmosphere (iNajita et al. 
2003h . 



Results 



3.1. Spitzer-IRS 



The IRS spectra of AS 205A and DR Tau show a large number of water emission lines 
(see Figured]), the vast majority pure rotational transitions with quantum numbers in the 
range 20-50. The i?=600 resolution of the SH module results in unresolved, and signifi- 
cantly blended, lines. Consequently, the identification with water is achiev ed via a simple, 



isoth ermal LTE ring model computed using the HITRAN 2004 database (iRothman et al. 



20051 ). The ring model consists of a single temperature (T) and column density (A^) gas 
in Keplerian rotation. Model parameters include T, A^, solid angle (fl), stellar mass (M^), 
inclination (i) and characteristic radius (r). In addition, the local line broadening is assumed 
to be Gaussian, with specified width (a). M^, and i were fixed to values taken from the litera- 
ture (j Clarke fc Bouvieiil2000l : iKitamura et al.ll2002l : lEisner et al.ll2005l : [Andrews fc Williams 
20071 ). The spectrally unresolved IRS data cannot completely constrain the remaining model 
parameters, but do require high T (T > 500 K) and A^ (A^ > 10^'' cm~^ for (T~2 km s^^). 
Figure [1] compares the IRS spectra with a model in which T and A^ are constrained using 
the NIRSPEC data described below, but where Q has been increased by a factor of 2. 

The addition of CO2 (and, for AS 205A, OH) improves our fits, as shown by the AS 
205A model in Fig. [H With T and Q fixed, we find H2O:CO2=40 for DR Tau and 20 for 
AS 205. Detections of these and other molecules will be discussed in detail elsewhere. 
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3.2. NIRSPEC 

While tropospheric water vapor blocks significant portions of incoming radiation near 
3 fim, there are a handful of rovibrational lines that are optically thin in the atmosphere but 
readily excited at temperatures typical of the inner ~ 1 AU of circumstellar disks. This is 
illustrated in Figure [H in which spectral regions with <70% atmospheric transmission have 
been removed, yet many prominent (line-to-continuum ~ 5- 10%) H2O emission features 
remain. These are primarily 2/3 = 1 ^ lines, with upper level energies of ~ 5,000 -10,000 
K. A few weaker (line-to-continuum ~ 1%) features are nominally detected from 3.13-3.18 
fim and are used as model constraints. We also detect the P4.5 (2.9344 /im; Tup ~ 5400 K) 
and P9.5 (3.1788 /xm; Tup ~ 7500 K) OH doublets. The M-band spectra include portions of 
the i^CO/^^CO v=1^0 and ^^CO v=2^1 P- and R-branches (Figure [3]). 

We fit all spectra with the ring models described in §3.1[ Because the CO lines span a 
large range of excitation energies, across which they transition from optically thick to opti- 
cally thin, N, T, and fl can be determined uniquely. Also, a constr aint on iVco is provide d 



by the ^CO lines (assuming a ^^CO. ^CO ratio, here fixed at 77, c.f. iBlake fc Boogertll2004j ). 



Thus, we began by fitting the CO emission b y comparison with an excitation diagram of the 



f=1^0 transitions (Fig. [3]), as described in ISalyk et al.l (120071 ). For our nominal H2O and 



OH fits, we then fixed Q and T to the CO-derived values to determine N - a reasonable, 
though not ideal, assumption, given the similarity in line shapes (see Fig. [3]). 

Our nominal ring model is shown in Figure [21 with best-fit parameters in Table [H and 
uncertainties of order AA^~30%, AT~15% and Af]~10%. We have fixed r = 3 AU to 
roughly match the emission profiles since the lineshape does not affect the integrated fiux; 
we discuss lineshapes below. Q corresponds to areas of 0.2 and 0.4 AU^ for DR Tau and AS 
205, respectively, at distances of 140 and 120 pc. Uncertainties in the nominal H2O and OH 
fits are linked to CO uncertainties, but if T is instead allowed to vary, the H2O spectra are 
consistent with T~900-1200/900-1100 K for DR Tau/AS 205A and ranges in of about 
1 order of magnitude. For OH, we can only constrain T ^^900 K. 

The model results depend crucially on the assumed local line width [a), which is un- 
known. For our nominal model, we have adopted a = 2 km s~^ — the sound speed for H2 at 
1000 K. However, we have also tested values between ~0.3 and 10 km s~^. We find that the 
opacity needs to remain similar from fit to fit, so that a change in requires a proportional 
change in N. To maintain the fiux level, the total number of molecules must remain similar, 
so an increase in a implies a decrease in Q. Accordingly, the C0:H20 ratio changes by no 
more than a factor of ~ 2 between models. (If ThjO is not fixed, however, this ratio may 
change). In addition, T is only weakly affected by a, since it is set by the overall shape of 
the spectrum, rather than by the absolute flux levels. 
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NIRSPEC spectrally resolves the molecular emission, allowing an investigation of the 
disk gas kinematics. In Fig. [3] we compare the unblended H2O transition at 2.931 fim with 
a CO line composite constructed from an average of all f =1^0 lines. In general, Keplerian 
disks should produce double-peaked emission profiles, with peak separation set by the inner 



and outer emission radii (nn. Tout), as well as the temperature gradient (IHorne fc Marsh 



19861 ). The CO and H2O line profiles, however, are more Gaussian, or even Lorentzian in 
shape, with relatively wide wings and a narrow peak. The narrow core could be produced by 
a ring with a very large r^^ (>10 AU) if the peak-to-peak separation is equal to the NIRSPEC 
resolution (given the and i in Table [1]). However, this is difficult to reconcile with the 
high T derived from the model fits and the observed line wings. 

Another way to produce a single peak would be to have the line emission extend to 
significant radii (rout of order 7 AU, for a temperature profile that is constant with radius), 
such that the contribution from low-velocity portions of the disk infills the center of the line 
profile. However, this model produces too much flux to match observations. Additionally, 
passive disk models predict a steep temperature decline with radius. A more promising 
solution may be that some flux is produced via fluorescence, in which molecules in the 
flared outer disk (typically out to ~5 AU) inte rcept and re-radiate inn er disk IR continuum. 



thus filling in the line profiles at low velocity (IBlake &: Boogertll2004l ). With canonical disk 
flaring and temperature structure, we find that a resonance fluorescence model with 100% 
scattering efficiency can indeed produce fluxes comparable to those of thermal emission; a 
more thorough investigation is left as future work. 

In any case, the amount of flux in the line profile wings is consistent with emission 
arising in a Keplerian disk with a small inner radius. With this assumption, the line wing 
velocity can be used to estimate ri^. Because turbulent motion and disk structure can both 
have non-negligible effects on the line shapes, we adopt here the velocity, fin, at 2.2 times 
the Half Width at Half Maximu m — a value foun d to be appropriate for the moderate- 



inclination disk around CM Aur (jSalyk et al.l 120071 ). To compute fin we use Gaussian fits 



to the CO and H2O lines discussed above, as well as the average of two-Gaussian fits to 
the OH emission doublets, deconvolved with the NIRSPEC instrument response function. If 
turbulent velocities are low, radii associated with fin will be overestimates, and so we also 
compute f at 3cr, which should represent a lower limit to Tin (see Tabled]). 



4. Discussion 



We have convincingly detected H2O emission in the 10-20 fim region with Spitzer-IRS 
as well as H2O and OH emission near 3 fim with NIRSPEC, arising from the disk atmospheres 
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surrounding two cTTs. By combining these data with constraints from spectrally resolved 
CO f=1^0 emission, we find that the excitation temperatures are typical of terrestrial 
planet forming regions (~ 1000 K). Additionally, line wing velocities imply inner emission 
radii no larger than 1 AU for all molecules. Therefore, further observations and analyses of 
these species may provide a new window into terrestrial planet forming regions. 

Observations of H2O in disks may also provide constraints on disk evolution and water 
transport. In a disk in which water vapor transport is controlled by diffusion, the presence 
of a condensation front ('snow line') can potentially dry out an inner disk in as little as 10^ 
yr, for typical disk viscosities (IStevenson &: Lunindll988l ). The snow line for the early solar 
syster a is estimated to be near 3 AU, but would be further out for these high accretion rate 
stars (ILecar et al.ll2006l ). Therefore, the presence of significant water vapor inside 1 AU may 
be evidence for inward radial migration, or upward mixing, of icy solids, followed by evap- 
oration. With both solid migration and diffusion controlling water transport, water vapor 
concentrations become sensitive to an array of parameters, including disk vis cosity and plan- 
etesimal growth rates, and water concentrations may trace disk evolution (jCiesla fc Cuzzi 
2006h . 



The detection of strong OH emission is also interesting, for OH is known to be an 
important ingredient to the chemistry in t he inner regions of pro toplanetary disks. It also 



acts as a coolant near the disk surface (e.g. iDuUemond et al.l 120071 ). and photolys is of OH to 



0(^p ) has been invoked to explain the observed strengths of the 6300 A O I line ( Acke et al. 
20051 ). Detections of OH in these and other disks will be discussed in greater detail in 
Mandell et al., in prep. 

An understanding of the local line broadening in disks will be crucial for obtaining 
accurate absolute column densities of molecular species, as smaller local line widths require 
lower column densities and vice versa. Nevertheless, we find that molecular ratios are more 
robust, and remain similar when utilizing different disk models with the nominal fits having 
C0:H20~ 10-15 and H20:OH~3-4, provided ThjO = ^oh = ^co- Stronger constraints on 
both T and will require high line-to-continuum observations across a wider wavelength 
range, but the available data reveal that H^Oi CO is similar to, or s lightly higher than that 
in the dense clouds out of which disks form (IBoogert et al.ll2004j ). Verti cally integrated 
H20: 0H ratios have been predicted to be much higher (> 10^) near 1 AU (IMarkwick et al. 



2OO2I ). but significa nt OH abundances are exp ected near the C^/C/CO transition zone at 
large scale heights (IKamp fc DuUemondl 120041 ) . Also, it is important to keep in mind that 



disk transport is not included in current chemical models. 



Our analysis of the L- and M-band spectra thus far has assumed thermal excitation 
of the vibrationally excited states, but this need not be the case. Af =1 emission can also 
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be driven by the stellar UV ou tput or by the absorpt i on of IR photons from the star or 
disk, followed by fluorescence (Idello Russo et al.l l2004l : iBlake &: BoogertI 120041 ). For OH, 
vibrationally excited emission caii also result from the (photo) dissociation of H2O followed 
by prompt emission (IBonev et al.ll2006l ). In fact, fluorescence may be necessary to explain 
observed line profile shapes (see §3.20 . and could explain the surprisingly strong OH emission. 
A better understanding of these processes may therefore prove essential for relating observed 
line fluxes to disk column densities. 

If emission in the line wings truly represents emission from the disk inner radius, rel- 
atively large inner emission radii are implied for CO, OH, and H2O. Both DR Tau and 
AS 205A have had their dust inner rim sizes measured with the Keck Interferometer. At 
0.11 ± 0.03 and 0.07 ± 0.01 AU , respectively, they are consistent with the location of dust 
sublimation ( lAkeson et al.ll2005l : Eisner et al.ll2005h . while the inner emission radii we derive 
are at least 3 times larger. Perhaps these high accretion-rate stars can stir up the inner disk 
enough to erase the vertical temperature gradients necessary for formation of emission lines. 
Another possibility may be that replenishment rates cannot keep up with photo-dissociation 
at the disk surface. 

The ease with which H2O and OH are now detected in the terrestrial planet-forming 
zones of certain disks holds the promise of more extensive and detailed studies. Although 
the ground based echelle observations presented here probe only the inner disk atmospheres, 
future high spectral resolution follow-up studies of Spitzer-detected water emission from 
protoplanetary disks will provide constraints at a variety of disk radii and vertical depths. 
Additionally, spectro-astrometric observation s are being used to direct ly constrain emitting 
location and structure of the molecular gas (jPontoppidan et al.ll2008l ). Eventually, a suite 
of constraints combined with non-LTE radiative transfer models will allow us to rigorously 
address such far-reaching questions as: How does the water vapor abundance vary within 
and between disks? What constraints can be put on models for the chemistry and transport 
of volatiles in disks, and what implications will these have for the early evolution of the Solar 
System? 
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14.0 14.5 15.0 15.5 16.0 16.5 17.0 

Wavelength [am] 

Fig. 1. — A portion of the Spitzer-IRS data (in black) and molecular emission models (in 
red). The three components of the AS 205A model are shown in color below. The largely 
featureless spectrum of CoKu Tau/3 is included to demonstrate the dynamic range achievable 
by the IRS and the level of systematics, such as residuals from the de-fringing process. For 
both spectra, splines have been used to remove broad continuum features from silicates. 
Line-to-continuum excesses are ~5-10%. 
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Fig. 2. — Comparison between portions of the NIRSPEC L-band data and best-fit H2O/OH 
disk models (offset). 
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Fig. 3.— Top, left: a portion of the NIRSPEC M-band spectra. ^^CO/^^CO v=1^0 tran- 
sitions are marked with sohd and dashed hnes, respectively. ^^CO v=2^1 transitions are 
marked with dot-dashed lines. Top, right: CO and H2O emission lines are overplotted. Bot- 
tom: an excitation diagram with la error bars and best-fit models to the v=l—>0 transitions. 
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Table 1: Fits to NIRSPEC data 





DR Tau 


AS 205A 


CO ^;in, V3a. [kms ^] 


27, 32 


35, 40 


COrm, ra. [AU] 


0.8, 0.6 


0.5, 0.4 


H2O ^;in, [kms"^] 


24, 28 


36, 42 


H2O rm, ra, [AU] 


1.0, 0.7 


0.4, 0.3 


OH Vin, [kms"^] 


28, 32 


34, 39 


OHn,, ra. [AU] 


0.7, 0.5 


0.5, 0.4 


r[AU] 


3 


3^ 


a [km s~^] ^ 


2 


2 


[Mq] 


0.76 


1.2 


i[°]- 


67 


47 


^co,H20,OH [sr] 


2 X 10-16 


6 X 10-16 


7cO,H20,OH [K] 


1000 


1000 


A^co [cm-2] 


8 X 10^^ 


9 X IQi^ 


A^HaO [cm-2] 


8 X 10^^ 


6 X 10^^ 


A^OH [cm-2] 


2 X 10^^ 


2 X 10^^ 



"Fixed 

''Chosen to approximately match hne shapes 



